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1 Problem 

I MOTIO MQTll 

This is %.j|^s^_p the series^ oj^t^ree papers (following pV, 18]) and the 
theorem li.i and corollary li.2 below constitute the final goal of this series. 



corollary li.2 below constitute the final goal 
Arguments of this paper are rather standard; all the heavy lifting was done 
before. Let us consider the following operator (quantum Hamiltonian) 



1<J<N l<j<k<N 



-1 



m 



(1.2) ^= ^ je, = ^^{M.\C^) 

l<n<N 

with 

(1.3) H"" = {{iV - A) ■ aY - V{x) 
Let us assume that 

(1.4) Operator H is self-adjoint on i^. 
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We will never discuss this assumption. We are interested in the ground state 
energy E^(/4) of our system i.e. in the lowest eigenvalue of the operator H 
on Sj: 

(1.5) E*„(0) = infSpecH on S3 

as A = and more generally in 

E;^ = inf ^inf SpeCj-^ H + i y \V x A\^ dx 



(1.6) 
where 



:i.7) 



X - x„ 



l<m<M 

A/^Z>1, Z:=Zi + ...Zm, 



Zi > 0, ... ,Zm > 



M is fixed, under assumption 



(1.9) 



< a < K*Z~^ 



thm-1-1 
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with sufficiently small constant n* > 0. 
Our purpose is to prove 

11-9 2 

Theorem 1.1. Under assumption (11.9) as N > Z — CZ 3 



(1.10) 

provided 
fl.ll) 



E; = 4'+ E 2Z^S(aZJ + 0(A/¥+aa-3A/2) 



l<m<M 



mm I 

l<m<m'<M 



| Lieb-l li vrii : ground 

where Sf/ is a Thomas-Fermi energy 4§ge fnTor I IS ) and S(Z„) 

, 1^ 



Z?, are 



cor-1-2 
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magnetic Scott correction terms (^see 

Combining with the properties of the Thomas-Fermi energy we arrive to 
Corollary 1.2. Let us consider x^r, = minimizing full energy 



(1.12) 

Assume that 
(1.13) 



E/v + ^ ] ZmZn 



-1 



l<m<m'<M 



z„ 



A/ 



Vm 



1-10, 



1 i-iy , 16 V 

Then a > N^^ and the remainder estimate in (11.10) is 0[N^) 
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Remark 1.3. As a = the remainder estimate fli.i2) was proven in ana 
the^remaipder estimate 0{^Ni[N~^ + in "FS] for atoms (M = 1) and 

\\2\ for M > 1; this better asymptotics contains also Dirac and Schwinger 
correction terms. Unfortunately I was not able to recover si^ch remainder 
estimate here unless a satisfies stronger assumption than (li.9). I still hope 
to achieve this better estimate without extra assumptions. 



p^^ satisfy equations 



Recall that Thomas-Fermi potential 1/1/^"^ and Thomas-Fermi density 



1^ (1.14) 
and 
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TF 



37r2 



(1/1/'^) 



:i.i5) 



1/1/ 



TF 



1 
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We prove theorem li.I m sections b "Lower estimate" and B „"Upper 

segt-C cor-1-2 

Estimate" . Section 4 "Miscellaneous" is devoted to corollary 11.2 and a brief 
discussion. 
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2 Lower estimate 

Consider corresponding to H quadratic form 
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(2.1) (Hvl/,vl/) = 5^(H°vl/,vl/) + ( |x,-x,|-^vl/,vl/) = 

j l<j<k<N 

5^(H,^vl/,vl/) + ((\/- l/l/)vl/,vl/) + ( ki-x/<r'M/,vl/) 

J l<j<k<N 



with 



2-2 



(2.2) 



H= ((/V-A)-o-)^- 1/1/ (x) 



where we select W later. By Lieb-Oxford inequality the last term is estimated 
from below: 



2-3] (2.3) 



l<j<k<N 
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where 



2^ (2.4) p^{x) = N j |vl/(x;x2,...,XA,)pc/x2---c/x/v 

is a spatial density associated with U' and 
2-^ (2.5) D(p,p') := |x-y|-V(x)p'(y)a'xc/y 



Therefore 
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(2.6) (Hvl/,vl/)> 

5^(H,^.vl/, VI/) - 2((\/ - l/l/)vl/,vl/) + D(p^, p,;)-C L|c/x = 

^(H^.vl/, VI/) - 2D(p, p^) + D(pvi., p^)-C I pI dx = 
J 

J2{H,^^, VI/) - D(p, p) + D(p - p^,p- p^) - C J pI dx 



as 

2^ (2.7) 1/1/- \/= |x|-Sp. 

Note that due to antisymmetry of v|/ 



^ (2.8) E(^><.^' ^) ^ E ^ Tr-(H) 

j l<j<N:\j<0 

where Xj are eigenvalues of H. _ 
T(i,^|timate the last term in (2.6) we reproduce the proof of Lemma YLS 
-33]: 



from 

According to magnetic Lieb-Thirring inequality for U > 0: 
2^ (2.9) ((^° - ^)^' ^) > - C / L/^/^ dx - C-f-^U"^ dx--f f B^dx 

j<N •' 

B = VxA, 7>0is arbitrary. Selecting U = (3 min(p^''~ \7pJ/^) with /3 > 
small but independent from 7 we ensure |l/p\|/ > CU^^^ + C'j~^U^ and then 



10] (2.10) / min(p5/^7p4/3■)^^_^ /" g2 
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which imphes 

(2.11) j pfdx < j min(pf , ^p'l^)dx + jj p^dx < 

c-f-^ ((H°)vl/,vl/) + c f B^dx + c-fN 

J:Xj<0 ' 

where we use / py^dx = N. 

| ES3 

Remark 2.1. As one can prove easily (see also l[ES3]) that 
(2.12) 



we conclude that 



2-13 (2.13) 



^((H°)vl/,vl/)<CZtA/ 

j<N 



pT^x <CZ-3N+Ci I B'dx 



It is sufficient unless we want to recover Dirac-Schwinger terms which 
unfortunately are too far away for us. 

Therefc^e skipping the non-negative third term in the right-hand expres- 
sion of (b.6) we conclude that 



2^(2.14) (Hv|/,v|/) + - /"iV X /^Pc/x > 

a J 

Jr{H) + {- - Q) f \V X A\^ dx - D{p, p) - CN'^ . 
a J 

Applying Theorem b.2 from |I8] we conclude that 



2-15 



(2.1^_)^|,he sum of the first and the second terms in the right-hand expression 
of (c2.14) is greater than 



15n'' 
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To prove this estimate one needs just to rescale x h-> xA/5, a i— a A/? and 
introduce /? = A/~,3,„ai|iji k = aN. Here one definitely needs the regularity 

^ut we have them as p = p^"^, 1/1/ = 1/1/^'^. Also one 
Ci" brings correction not exceeding C2aN^ as aN < 1. 
Meanwhile for p = p^^ , W = W'^^ 



properties like in 
can see easily that 



(2.16) 



J W'2 dx - D{p, p) = . 



I thm-l-l 

Lower estimate of Theorem li.i has been proven. 

Remark 2.2. P,^^^^^ ■, ^ = W^^ delivers th^_]L)iaximum of the right-hand 
expression of (I2.i6) among p, 1/1/ satisfying ([2.7). 
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3 Upper Estimate 

li vrii : p:j:ound 

Upper estimate is easy. Plugging as the Slater determinant (see IS f.e. 

oi Tpi, ... ,%Ijn where ipi, ... ,ipf\i are eigenfunctions of Ha,w we get 



3:^(3.1) (Hvl/,vl/) =Tr(H4,^v-AA/) + AA/A/+ 

y"(l/l/- V){x)p^i,{x) dx + D{p^, p^y 



^N{N-1) II |xi-X2|-Xxi,X2;x3 x^)\^ dx^ ■ ■ ■ dx^ 



2 

where we don't care about last term as we drop it (again because we cannot 
get sharp enough estimate) and the first term in the second line is in fact 

3^ (3.2) -2D(p,p^); 

provided holds. Thus we get 



3-3 



(3.3) Tr {Ha,w-Xn) + XnN-D{p,p) + D{pm,-p,Ph,-p)+ 

1 



IdAl^dx 



where we added magnetic energy. Definitely we have several problems here: 
A/v depends on A and there may be less than N negative eigenvalues. 
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However in the latter case we can obviously replace N by the lesser 
number N := max(n < ,^_^^n < 0) ss E ^ is decreasing function of A/. In this 
case the first term in (S.3) "woyM.^^'S just Tr~{HA,w) and the second would 
be 0. Then we apply theory of IjlS] immediately without extra complications. 

Consider A a minimizer (or its mollification) for potential W = W^^ 

2 

and < 0. Then with an error 0(A/3) 

(3.4) #{A, </i} = J{W-fi)ldx+0{Nl). 

ne c^p prove easily using the regularity properties of A establishe|l 
L |lS] and the same rescaling as before. Note that the first term in (b.4) 



O] 
in 



differs from the same expression with fi = (which is equal to Z) by 
X /i(A/^/^)^/^ ■ = jjiN'^/^. Then as the left-hand expression equals A/, 
and A/ - Z = 0(A/i, we conclude that |A/v| = 0{N). 

There|o^e modulo 0{^N^ + fi:a~^A/^) the sum of the first and the last 
term in (B.3) is equal to 

(3.5) il^ /( - ^'^)+ + E 2^^5(«:Z^) 



and modulo 0{N—^\l) = 0{NI) one can rewrite the first term here as 

(3.6) i^Jmdx-x.^, 

and "^^t^ the same error the second term here cancels term A a/ A/ in 



then (3.3) becomes 



(3.7) ^ [^ldx + J2 2Z^S(«:Z„) - D(p, p) + D{p^ - p, - p) 

and as 1/1/ = l/l/""""^, p = p^^ the first and the third term together are £'^^ , so 
we get again 8'^^ + Y.m 2Z^S{KZm). ^ 
Now we need to estimate properly the last term in (cB.V) i.e. 

(3.8) lJJ\x- yr'iP^x) - P'^'ix)) (P4y) - P'^'iy)) dxdy. 

1 1-14 

Rescaling as before, and using (11.14) we conclude that it does not exceed 



,5 



(3.9) A/3 / / g{xyg{yyr'{x)r'{y)\x - y\-' dxdy 



where ^ is p of I lij and we know that g = as i < 1 and g = as £ > 1. 

Estimating integral by the (double) sum of integral as i{x) < 1, i{y) < 1 
and £(x) > 1, £(y) > 1 we get (increasing C) 



and 



C I \x — y\ "^|x| ^|y| ^ c/yc/x X 1 

'{|y|<W<i} 



c k~y| "^1^1 ^|y| ^ c/yc/x X 1 

J{\y\>\x\>i} 



respectively. i thm-l-l 

This concludes the proof of the upper estimate in Theorem II. 1 which is 
proven now. 
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4 Miscellaneous 

| cor-l-2 

Proof. Proof of corollary 11.2 Optimization with respect to xi, 



, X/V7 implies 



(4.1) 



l<m<M 



where E* = E*(xi, ... ,xm; Zi, ... , Z„, N) an ^^_^_= ^*i^m, Zm) are calculated 
for separate atoms. In virtue of theorem li . 1 aiid (li.9) then 



(4.2) 



l<m<M 



however due to strong non-binding theorem in Thomas-Fermi theory (see f.e. 
I Sol ^ ^_ _i | 4-g ^ 

l[S]) the left-hand expression is x a as a > N 3 and therefore (kl.2) implies 



"4^ (4.3) 



rem-4-1 



a > eiA/ s 



and a"^A/ < A/f . 

On the other hand, there is no binding with a < A/~3 because remainder 
estimate is (better than) CN^ and binding energy excess is x A/3. □ 

Remark 4-1- Similar arguments work if we improve A/t to N'^ with > | 
but without improving a~^N part of the remainder estimate we would not 
pass beyond 0{N^). 
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There are several questions which after I[I8] could be answered in this 
framework by the standard arguments with certain SjJPf^but we postpone 
it, hoping to improve remainder estimate 0(/?~3) in : 



problem-4-2 Problem 4-2. (i) Investigate case N < Z — CZi\ 

(ii) Estimate from above excess negative charge (how many extra electrons 
can and bind atom) ionization energy (E^_]^ — E^); 

(iii) Estimate from above excess positive charge in the case of binding of 
several atoms i.e. estimate Z — A/ as 



4-4 



(4.4) E*(xi, 



xm;Zi, ... ,Z„, A/) < min E;;(Z^, A/^). 

A/l I^m- 

Ni + ...+Nm^N 
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